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Atomic physics and hadron physics are both based on Yang Mills gauge theory; in fact, quantum 
electrodynamics can be regarded as the zero-color limit of quantum chromodynamics. I review a 
number of areas where the techniques of atomic physics provide important insight into the theory 
of hadrons in QCD. For example, the Dirac-Coulomb equation, which predicts the spectroscopy 
and structure of hydrogenic atoms, has an analog in hadron physics in the form of light-front 
relativistic equations of motion which give a remarkable first approximation to the spectroscopy, 
dynamics, and structure of light hadrons. The renormalization scale for the running coupling, 
which is unambiguously set in QED, leads to a method for setting the renormalization scale in 
QCD. The production of atoms in flight provides a method for computing the formation of hadrons 
at the amplitude level. Conversely, many techniques which have been developed for hadron physics, 
such as scaling laws, evolution equations, and light-front quantization have equal utility for atomic 
physics, especially in the relativistic domain. I also present a new perspective for understanding the 
contributions to the cosmological constant from QED and QCD. 

I. INTRODUCTION 

Quantum Electrodynamics, the fundamental theory of leptons and photons which underlies all of atomic and 
molecular physics, and Quantum Chromodynamics, the quark and gluon theory with three colors underlying hadronic 
and nuclear physics, are both derived from Yang-Mills gauge theory. The Yang-Mills Lagrangian for SU{Nc) is 
invariant under arbitrary color rotations and phases at each point of space and time. In fact, in the limit where 
the number of colors Nc vanishes, with asCp = held fixed {Cp = {N^ — l)/2Nc)^ QCD becomes equivalent to 
Abelian gauge theory [1]. This analytic connection as a function of Nc between QCD and QED provides a valuable 
link between the two fields; processes and analyses in QCD must connect at zero color to the analogous reactions and 
procedures of QED. 

In this paper, I will review a number of areas where the techniques of atomic physics give important insight into 
the theory of hadrons, the color-singlet bound states of quarks and gluons in QCD. For example, the Dirac-Coulomb 
equation, which predicts the spectroscopy and structure of hydrogenic atoms has an analog in hadron physics in 
the form of relativistic frame-independent equations of motion derived from light-front holography [2] which give a 
remarkable first approximation to the spectroscopy, dynamics, and structure of light hadrons. The renormalization 
scale for the running coupling which is unambiguously set in QED leads to a solution for setting the renormalization 
scale in QCD. The production of atoms in flight provides a method for computing the formation of hadrons at the 
amplitude level. Conversely, many techniques and theorems developed for hadron physics, such as scaling laws, 
evolution equations, and light-front quantization have equal utility for atomic physics, especially in the relativistic 
domain. 



II. PRODUCTION OF EXOTIC ATOMS IN FLIGHT AND HADRONIZATION AT THE AMPLITUDE 

LEVEL 

Relativistic antihydrogen was first produced in 1995 at CERN-LEAR [3] and at the Fermilab Antiproton Accu- 
mulator [4]. The production mechanism [5] is illustrated in fig. 1 (a). The incident antiproton beam produces a 
Bethe-Heitler electron-positron pair in the Coulomb field of a target nucleus pZ pe^e~ Z [pe~^]Z. The comoving 
off-shell p and e+ then coalesece into antihydrogen atoms via the Schrodinger Coulomb wavefunction which connects 
the off-shell state to the on-shell anti-atom. The atom is dominantly in its IS ground state. In principle, one can 
measure its "anti-Lamb- Shift" using the Robiscoe level-crossing method [6]. 

The production of antihydrogen in flight provides important insight into the dynamics of hadron production in 
QCD. For example, the A{sud) baryon can be produced at high longitudinal momentum fraction xp in pp ^ AX 
reactions by the coalescence of the ud valence quarks of the beam with a comoving strangeness quark. This method 
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can be generalized to produce heavy hadrons such as Kc{cud)^ Kij^ double charmed baryons, etc., using the high x 
intrinsic heavy quarks which exist in the higher Fock states of the proton wavefunction [7]. 

The analog of intrinsic charm in hadrons is the content of positronium. The |e+e~/i+/i~ > Fock state appears 

through the cut of the muon-loop light-by-light contribution to the self energy of the positronium eigenstate. In this 
Fock state, the muons carry almost all of the momentum of the moving atom since the off-shell virtuality is minimal 
at equal velocity. In QED the probability for intrinsic leptons LL exist in positronium scales as ^/m\ whereas in 
QCD the probability of intrinsic heavy quarks in the wavefunction of a light hadron scales as 1/mg because of its 
non-Abelian couplings [8, 9]. 

The production of a qq meson in an e^e~ annihilation event is illustrated in fig. 1(b). One first calculates the 
T matrix element for the production of off-shell quarks and gluons at the amplitude level using light-front time- 
ordered perturbation theory. The light-front wavefunction of the meson then converts the off-shell comoving qq pair 
into the final-state meson. The confined colored quarks thus never appear on-shell. This first-principle method for 
forming hadrons in QCD [10] can replace phenomenological jet hadronization models such as PYTHIA. The light- 
front wavefunction required for calculating "hadronization at the amplitude level" [10, 11] is the frame-independent 
analog of the Schrodinger wavefunction of atomic physics. It is obtained from the eigensolution of the QCD light- 
front Hamiltonian quantized at fixed light-front time r which can be determined by solving the Heisenberg matrix 

using a method such as discretized light-cone quantization (DLCQ) [12] or using the 
AdS/QCD approach together with Light-Front Holography [2]. 

It is very interesting to produce "true muonium", the [/i+/i~] bound state Lebed and I [13] have discussed QED 
production and decay mechanisms, such as electroproduction of relativistic true muonium below the threshold 
via e~Z [/i+/i~]e~Z or e+e~ [/i+/i~]7. See fig. 2. The APEX electroproduction experiment [14], which will 
search for dark matter candidates at Jefferson Laboratory, could be the first to see this exotic atom. Studying the 
precision spectroscopy of the [/i+/i~] atom is important in view of the anomalies seen in the muon g — 2 [15] and the 
fi~p Lamb shift [16]. 

"Atomic Alchemy "refers to the transition between a muonic atom into an electronic atom: {/J^~ Z) {e~ Z)uei^^ 
via the weak decay of the bound muon and the subsequent capture of its decay electron. Greub, Wyler, Munger 
and I [17] have shown that such processes provide a laboratory for studying the relativistic high momentum tail of 
wavefunctions in atomic physics; in addition, they provide a simple toy model for investigating analogous exclusive 
heavy hadronic decays in quantum chromodynamics such as 5 ^ ireu. 

The QCD analog of a molecule in QCD is a bound state of heavy quarkonium with a nucleus such as [J/tJjA] [18, 19]. 
The binding occurs through two-gluon exchange, the hadronic analog of the Van der Waals interaction. Since the 
kinetic energy of the J/i/j and the nucleus are both small, one expects to find produce these exotic hybrid states 
at threshold. Examples of nuclear-bound quarkonium are the \uuduudss > and \uuduudcc > resonances which 
apparently appear as intermediate states in pp pp elastic exchange. These resonances can account [20] for the large 
spin-spin Ann correlations [21] observed at the strangeness Ecm — 3 GeV and Ecm — 5 GeV and charm thresholds. 

At high energies, Compton scattering on an atom ^yA — > ^yA is dominated by the Thomson amplitude - the elastic 
scattering of the photon on the atomic electrons. The analog in hadron physics is the scattering of photons on quarks 
7Q ~^ via a local seagull or instantaneous light-front term or which gives an energy-independent contribution to 
the Compton amplitude proportional to the charge squared of the struck quark - a contribution which has no analog 
in hadron scattering reactions. Llanes-Estrada, Szczepaniak, and I [22] have shown that this local contribution has 
a real phase and is universal, giving the same contribution for real or virtual Compton scattering for any photon 
virtuality and skewness at fixed momentum transfer squared t. The t-dependence of this J = fixed Regge pole is 
parameterized by a yet unmeasured even charge-conjugation form factor of the target nucleon. The t = limit gives 
an important constraint on the dependence of the nucleon mass on the quark mass through the Weisberger relation. 
The same J = amplitude enters the two-photon exchange contribution to muon-proton scattering, and thus also 
could contribute an important contribution to the jiip Lamb Shift. 

III. RENORMALIZATION SCALE SETTING 

A key difficulty in making precise perturbative predictions for QCD is the uncertainty in determining the renormal- 
ization scale /i of the running coupling Qfs(/i^). In the standard Gell-Mann-Low scheme for QED, the renormalization 
scale is simply the virtuality of the virtual photon [23]. Although the initial choice of renormalization scale to is 
arbitrary, the final scale t which sums the vacuum polarization corrections is unique and unambiguous. The resulting 
perturbative series is identical to the conformal series with zero /3- function. In the case of muonic atoms, the mod- 
ified muon-nucleus Coulomb potential is precisely —Za{—q'^)/q'^; i.e., /i^ = —q^. Again, the renormalization scale 
is unique. The same principle underlying renormalization scale-setting in QED for Nc = must also hold in QCD 
since the np terms in the QCD (3 function have the same role as the lepton Ni vacuum polarization contributions in 
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FIG. 1: (a) Production of Relativistic Antihydrogen. (b) Hadronization at the Amplitude Level in Electron-Positron Annihi- 
lation 
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FIG. 2: Production of True Muonium 

QED. Thus the same scale-setting procedure must be apphcable to ah renormahzable gauge theories. 

The purpose of the running couphng in any gauge theory is to sum ah terms involving the /3 function; in fact, 
when the renormalization scale \i is set properly, all non-conformal /3 7^ terms in a perturbative expansion arising 
from renormalization are summed into the running coupling. The remaining terms in the perturbative series are 
then identical to that of a conformal theory; i.e., the theory with /3 = 0. The divergent "renormalon" series of order 
OL^P^n\ does not appear in the conformal series. Thus as in quantum electrodynamics, the renormalization scale \i is 
determined unambiguously by the "Principle of Maximal Conformality (PMC)" [24, 25]. This is also the principle 
underlying BLM scale setting [26] An important feature of the PMC is that its QCD predictions are independent of 
the choice of renormalization scheme. The PMC procedure also agrees with QED scale-setting in the Nq limit. 



IV. LIGHT-FRONT QUANTIZATION 



The distributions of electrons within an atom are determined in QED using the Schrodinger wavefunction, the 
eigenfunction of the QED Hamiltonian. In principle, one could calculate hadronic spectroscopy and wavefunctions by 
solving for the eigenstates of the QCD Hamiltonian: = E\^) at fixed time t. However, this traditional method - 

called the "instant" form" by Dirac, [27] is plagued by complex vacuum and relativistic effects, as well as by the fact 
that the boost of such fixed-t wavefunctions away from the hadron's rest frame is an intractable dynamical problem. 
However, there is an extraordinarily powerful non-perturbative alternative - quantization at fixed light-front (LF) 
time r = t + z/c = x+ =x^ + x^- the "front-form" of Dirac. [27] In this framework each hadron R is identified as an 
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eigenstate of the QCD Hamiltonian H^f^I'^h) = Mfjl^n). where H^^^ = P^P^ = P'P^-Pl is derived directly 
from the QCD Lagrangian or action. The eigenvalues of this Heisenberg equation give the complete mass spectrum of 
hadrons. The eigensolution 1^//) projected on the free Fock basis provides the set of valence and non- valence light- 
front Fock state wavefunctions ^^///(x^, /c^i, A^), which describe the hadron's momentum and spin distributions and 
the direct measures of its structure at the quark and gluon level. If one quantizes the gluon field in light-cone gauge 
^4+ = -\- = 0, the gluons have physical polarization = ±1, there are no ghosts, so that one has a physical 
interpretation of the quark and gluon constituents. The constituents of a bound state in a light-front wavefunction 
are measured at the same light-front time r - along the front of a light- wave, as in a flash picture. In contrast, the 
constituents of a bound state in an instant form wavefunction must be measured at the same instant time t - - this 
requires the exact synchrony in time of many simultaneous probes. 

A remarkable feature of LFWFs is the fact that they are frame independent; i.e., the form of the LFWF is 
independent of the hadron's total momentum P+ = P^ -\- P^ and P±. The boost invariance of LFWFs contrasts 
dramatically with the complexity of boosting the wavefunctions deflned at flxed time t. [28] Light-front quantization 
is thus the ideal framework to describe the structure of hadrons in terms of their quark and gluon degrees of freedom. 
The constituent spin and orbital angular momentum properties of the hadrons are also encoded in the LFWFs. The 
total angular momentum projection [29] = X]r=i '^i '^'^7=1 conserved Fock-state by Fock-state and by every 

interaction in the LF Hamiltonian. The constituent spin and orbital angular momentum properties of the hadrons 
are thus encoded in their LFWFs. The empirical observation that quarks carry only a small fraction of the nucleon 
angular momentum highlights the importance of quark orbital angular momentum. In fact the nucleon anomalous 
moment and the Pauli form factor are zero unless the quarks carry nonzero . 

Hadron observables, e.g., hadronic structure functions, form factors, distribution amplitudes, GPDs, TMDs, and 
Wigner distributions can be computed as simple convolutions of light-front wavefunctions (LFWFs). For example, one 
can calculate the electromagnetic and gravitational form factors < p-\- q\j^{0)\p > and < p -\- q\t^^ {0)\p > of a hadron 
from the Drell-Yan-West formula - i.e., the overlap of LFWFs. The anomalous gravitomagnetic moment B{0) deflned 
from the spin-flip matrix element < p -\- q\t^^{0)\p > at g ^ vanishes - consistent with the equivalence theorem 
of gravity. In contrast, in the instant form, the overlap of instant time wavefunctions is not suflicient. One must 
also couple the photon probe to currents arising spontaneously from the vacuum which are connected to the hadron's 
constituents. The Light-Front method is directly applicable for describing atomic bound states in both the relativistic 
and nonrelativistic domains; it is particularly useful for atoms in flight since the LFWFs are frame-independent. It 
also satisfles theorems such as cluster decomposition. 

One can solve the LF Hamiltonian problem for theories in one-space and one-time by Heisenberg matrix diagonal- 
ization. For example, the complete set of discrete and continuum eigensolutions of mesons and baryons in QCD(1+1) 
can be obtained to any desired precision for general color, multiple flavors, and general quark masses using the dis- 
cretized light-cone quantized (DLCQ) method. [30, 31] The DLCQ approach can in principle be applied to QED(3+1) 
and QCD(3+1); however, in practice, the huge matrix diagonalization problem is computational challenging. 

V. ADS/QCD LIGHT-FRONT HOLOGRAPHY 

Recently a new nonperturbative QCD approach has been developed which leads to an elegant analytical and 
phenomenologically compelling flrst approximation to the full LF Hamiltonian method- "Light-Front Holography" . [2] 
Light front holographic methods allow one to project the functional dependence of the wavefunction ^{z) computed 
in the AdS flfth dimension to the hadronic frame-independent light-front wavefunction il^{xi^h^i) in 3 + 1 physical 
space-time. The variable z maps to a transverse LF variable C(^i7^±i)- The result is a single- variable light-front 
Schrodinger equation which determines the eigenspectrum and the LFWFs of hadrons for general spin and orbital 
angular momentum. The transverse coordinate is closely related to the invariant mass squared of the constituents 
in the LFWF and its off-shellness in the LF kinetic energy, and it is thus the natural variable to characterize the 
hadronic wavefunction. In fact is the only variable to appear in the relativistic light-front Schrodinger equations 
predicted from holographic QCD in the limit of zero quark masses. The coordinate z in AdS space is thus uniquely 
identifled with a Lorentz-invariant coordinate ^ which measures the separation of the constituents within a hadron at 
equal light-front time. 

The hadron eigenstates generally have components with different orbital angular momentum; e.g., the proton 
eigenstate in LF holographic QCD with massless quarks has L = and L = 1 light-front Fock components with equal 
probability. Higher Fock states with extra quark-anti quark pairs also arise. The resulting LFWFs then lead to a 
new range of hadron phenomenology, including the possibility to compute the hadronization of quark and gluon jets 
at the amplitude level. The soft-wall model also predicts the form of the non-perturbative effective coupling and its 
/3-function. [32] 



6 



Single/- ^On/ 
ouyyvnmetry 



' Sp-qxpq 




Leading-TVist 
Sivers Effect 



proton 



QCD S- and P- 
Coulomh Phases 
--Wilson Line 

''Lensing Effect" 



s>Iq\a. reversal 



FIG. 3: Origin of the Sivers single-spin asymmetry in deep inelastic lepton scattering. 



VI. LENSING AND THE SIVERS EFFECT 



A well-known phenomenon in QED rescattering via final-state Coulomb interactions. Although the Coulomb phase 
for a given partial wave is infinite, the interference of Coulomb phases arising from different partial waves leads to 
observable effects. 

The calculation of the Sivers single-spin asymmetry in deep inelastic lepton scattering in QCD is illustrated in fig. 3. 
The analysis requires two different orbital angular momentum components: 5- wave with the quark-spin parallel to the 
proton spin and P-wave for the quark with anti-parallel spin; the difference between the final-state "Coulomb" phases 
leads to di S ■ qx p correlation of the proton's spin with the virtual photon-to-quark production plane [33]. Thus, as 
it is clear from its QED analog, the final-state gluonic interactions of the scattered quark lead to a T-odd non-zero 
spin correlation of the plane of the lepton-quark scattering plane with the polarization of the target proton [33] . This 
leading-twist Bjorken-scaling "Sivers effect" is nonuniversal since QCD predicts an opposite-sign correlation [34, 35] 
in Drell-Yan reactions due to the initial-state interactions of the annihilating antiquark. The S— and P-wave proton 
wavefunctions also appear in the calculation of the Pauli form factor quark-by-quark. Thus one can correlate the 
Sivers asymmetry for each struck quark with the anomalous magnetic moment of the proton carried by that quark [36], 
leading to the prediction that the Sivers effect is larger for positive pions. The physics of the lensing dynamics involves 
nonperturbative quark-quark interactions at small momentum transfer, not the hard scale of the virtuality of the 
photon. It would interesting to see if the strength of the soft initial- or final- state scattering can be predicted using 
the confining potential of AdS/QCD. 
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VII. VACUUM CONDENSATES AND THE COSMOLOGICAL CONSTANT 

It is important to distinguish two very different concepts of the vacuum in quantum field theories such as QED 
and QCD. The vacuum is normally defined as the lowest energy eigenstate of the instant-form Hamiltonian - the 
vacuum defined by quantizing at fixed time t. In QED, the instant-time vacuum is saturated with quantum loops of 
leptons and photons. In calculations of physical processes one must normal-order the vacuum and divide the ^'-matrix 
elements by the disconnected vacuum loops. In contrast, the front-form (light-front) vacuum is defined as the lowest 
mass eigenstate of light-front Hamiltonian defined by quantizing at fixed r = t — z/c. The vacuum is remarkably 
simple in light-front quantization because of the restriction /c+ > 0. For example QED vacuum graphs such as e+e~7 
do not arise. The LF vacuum thus coincides with the vacuum of the free LF Hamiltonian. The front-form vacuum 
and its eigenstates are Lorentz invariant; whereas the instant form vacuum depends on the observer's Lorentz frame. 
The instant-form vacuum is a state defined at the same time t at all spatial points in the universe. In contrast, 
the front-from vacuum only senses phenomena which are causally connected; i.e., or within the observer's light-cone. 
Causality in quantum field theory follows the fact that commutators vanish outside the light-cone. In fact in the LF 
analysis the spatial support of QCD condensates is restricted to the interior of hadrons, physics which arises due to the 
interactions of confined quarks and gluons. The condensate physics is replaced by the dynamics of higher non- valence 
Fock states as shown by Casher and Susskind. [37] In particular, chiral symmetry is broken in a limited domain 
of size I/ttIt^, in analogy to the limited physical extent of superconductor phases. This novel description of chiral 
symmetry breaking in terms of "in-hadron condensates" has also been observed in Bethe-Salpeter studies [38, 39]. The 
usual argument for a quark vacuum condensate is the Gell-Mann-Oakes-Renner formula: = —2mq{0\qq\0)/f^. 
However, in the Bethe-Salpeter and light-front formalisms, where the pion is a qq bound-state, the GMOR relation is 
replaced by = —2mq{0\qjsq\7r)/fT^^ where p^^ = — (0|^75g'|7r) represents a pion decay constant via an an elementary 
pseudoscalar current. 

The cosmological constant measures the matrix element of the energy momentum tensor T^^ in the background 
universe. It corresponds to the measurement of the gravitational interactions of a probe of finite mass; it only senses 
the causally connected domain within the light-cone of the observer. If the universe is empty, the appropriate vacuum 
state is thus the LF vacuum since it is causal. One automatically obtains a vanishing cosmological constant from the 
LF vacuum. Thus, as argued in Refs. [40-42] the 45 orders of magnitude conflict of QCD with the observed value 
of the cosmological condensate is removed, and a new perspective on the nature of quark and gluon condensates in 
QCD is thus obtained. [40-42]. 

Acknowledgments 

I am grateful to the organizers of EXA2011 for their invitation to this outstanding interdisciplinary conference. I 
thank all of my collaborators whose work has been cited in this report, particularly Guy de Teramond, Rich Lebed, 
Ivan Schmidt, and Dae Sung Hwang. This research was supported by the Department of Energy contract DE-AC02- 
76SF00515. 



[1] S. J. Brodsky and P. Huet, Phys. Lett. B 417, 145 (1998) [arXiv:hep-ph/9707543]. 

[2] G. F. de Teramond and S. J. Brodsky, Phys. Rev. Lett. 102, 081601 (2009) [arXiv:0809.4899 [hep-ph]]. 

[3] G. Baur et al, Phys. Lett. B 368, 251 (1996). 

[4] G. Blanford, D. C. Christian, K. Gollwitzer, M. Mandelkern, C. T. Hunger, J. Schultz and G. Zioulas [E862 Collaboration], 

Phys. Rev. Lett. 80, 3037 (1998). 

[5] C. T. Hunger, S. J. Brodsky and L Schmidt, Phys. Rev. D 49, 3228 (1994). 

[6] R. T. Robiscoe, Phys. Rev. 138, A22 (1965). 

[7] S. J. Brodsky, C. Peterson and N. Sakai, Phys. Rev. D 23, 2745 (1981). 

[8] S. J. Brodsky, J. C. Collins, S. D. Ellis, J. F. Gunion and A. H. Hueller, 

[9] H. Franz, H. V. Polyakov and K. Goeke, Phys. Rev. D 62, 074024 (2000) [arXiv:hep-ph/0002240]. 

[10] S. J. Brodsky, G. de Teramond and R. Shrock, MP Conf. Proc. 1056, 3 (2008) [arXiv:0807.2484 [hep-ph]]. 

[11] S. J. Brodsky and G. F. de Teramond, arXiv:0901.0770 [hep-ph]. 

[12] H. C. PauU and S. J. Brodsky, Phys. Rev. D 32, 1993 (1985). 

[13] S. J. Brodsky and R. F. Lebed, Phys. Rev. Lett. 102, 213401 (2009) [arXiv:0904.2225 [hep-ph]]. 

[14] S. Abrahamyan et al [APEX Collaboration], Phys. Rev. Lett. 107, 191804 (2011) [arXiv:1108.2750 [hep-ex]]. 

[15] G. W. Bennett et al [Huon G-2 Collaboration], Phys. Rev. D 73, 072003 (2006) [arXiv:hep-ex/0602035]. 



8 



[16] R. Pohl, A. Antognini, F. Nez, F. D. Amaro, F. Biraben, J. M. R. Cardoso, D. S. Covita and A. Dax et al, Nature 466, 
213 (2010). 

[17] C. Greub, D. Wyler, S. J. Brodsky and C. T. Hunger, Phys. Rev. D 52, 4028 (1995) [arXiv:hep-ph/9405230]. 

[18] S. J. Brodsky, L A. Schmidt and G. F. de Teramond, Phys. Rev. Lett. 64, 1011 (1990). 

[19] M. E. Luke, A. V. Manohar and M. J. Savage, Phys. Lett. B 288, 355 (1992) [arXiv:hep-ph/9204219]. 

[20] S. J. Brodsky and G. F. de Teramond, Phys. Rev. Lett. 60, 1924 (1988). 

[21] G. R. Court et al, Phys. Rev. Lett. 57, 507 (1986). 

[22] S. J. Brodsky, F. J. Llanes-Estrada and A. P. Szczepaniak, Phys. Rev. D 79, 033012 (2009) [arXiv:0812.0395 [hep-ph]]. 
[23] M. Geh-Mann and F. E. Low, Phys. Rev. 95, 1300 (1954). 
[24] S. J. Brodsky and L. Di Giustino, arXiv: 1107.0338 [hep-ph]. 
[25] S. J. Brodsky and X. G. Wu, arXiv:1111.6175 [hep-ph]. 

[26] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983). 

[27] P. A. M. Dirac, Rev. Mod. Phys. 21, 392 (1949). 

[28] S. J. Brodsky and J. R. Primack, Annals Phys. 52, 315 (1969). 

[29] S. J. Brodsky, D. S. Hwang, B. Q. Ma and L Schmidt, Nucl. Phys. B 593, 311 (2001) [arXiv:hep-th/0003082]. 

[30] H. C. PauU and S. J. Brodsky, Phys. Rev. D 32, 2001 (1985). 

[31] K. Hornbostel, S. J. Brodsky and H. C. PauH, Phys. Rev. D 41, 3814 (1990). 

[32] S. J. Brodsky, G. F. de Teramond and A. Deur, Phys. Rev. D 81, 096010 (2010) [arXiv: 1002.3948 [hep-ph]]. 
[33] S. J. Brodsky, D. S. Hwang and L Schmidt, Phys. Lett. B 530, 99 (2002) [arXiv:hep-ph/0201296]. 
[34] J. C. Cohins, Phys. Lett. B 536, 43 (2002) [arXiv:hep-ph/0204004]. 

[35] S. J. Brodsky, D. S. Hwang and I. Schmidt, Nucl. Phys. B 642, 344 (2002) [arXiv:hep-ph/0206259]. 
[36] Z. Lu and L Schmidt, Phys. Rev. D 75, 073008 (2007) [arXiv:hep-ph/0611158]. 
[37] A. Casher and L. Susskind, Phys. Rev. D 9, 436 (1974). 

[38] P. Maris, C. D. Roberts and P. C. Tandy, Phys. Lett. B 420, 267 (1998) [arXiv:nucl-th/9707003]. 

[39] P. Maris and C. D. Roberts, Phys. Rev. C 56, 3369 (1997) [arXiv:nucl-th/9708029]. 

[40] S. J. Brodsky and R. Shrock, Phys. Lett. B 666, 95 (2008) [arXiv:0806.1535 [hep-th]]. 

[41] S. J. Brodsky and R. Shrock, Proc. Nat. Acad. Sci. 108, 45 (2011) [arXiv:0803.2554 [hep-th]]. 

[42] S. J. Brodsky and R. Shrock, Proc. Nat. Acad. Sci. 108, 45 (2011) [arXiv:0905.1151 [hep-th]]. 



